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Abstract 

We study the stability of the kaon condensed color-flavor locked (CFL) phase of dense quark 
matter with regard to the formation of a non-zero Goldstone boson current. In the kaon condensed 
phase there is an electrically charged fermion which becomes gapless near [is ~ 1.35 A and a 
neutral fermion which becomes gapless near fi^ ~ 1.61 A. Here, n s = m1/{2pp) is the shift in the 
Fermi energy due to the strange quark mass m s and A is the gap in the chiral limit. The transition 
to the gapless phase is continuous at [i s and first order at . We find that the magnetic 

(2) 

screening masses are real in the regime /i s < fi s , but some screening masses are imaginary for 
fi s > [i'P. We show that there is a very weak current instability for fi s > /xi 1 ^ and a more robust 

(2) 

instability in a small window near fj, s . We show that in the Goldstone boson current phase all 
components of the magnetic screening mass are real. There is a range of values of /i s below 2A in 
which the magnetic gluon screening masses are imaginary but the phase is stable with respect to 
electrically neutral fluctuations of the gauge field. 
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I. INTRODUCTION 



Understanding the response of a paired fermion system to an external pair breaking field 
has been an area of intense theoretical and experimental research over the last couple of 
years. In dense quark matter the problem arises in connection with the effect of a non- 
zero strange quark mass, see [l] for recent reviews. If quarks are massless and the flavor 
symmetry is exact then the ground state of three flavor quark matter at high baryon density 
is the color-flavor-locked (CFL) phase . The CFL state involves pairing between up 

and down, up and strange, as well as down and strange quarks. In the real world the up and 
down quarks are approximately massless but the mass of the strange quark is comparable 
to the QCD scale parameter. At high baryon density the effect of the strange quark mass is 
governed by the shift \x q = m 2 /{2^i) of the Fermi energy due to the mass. 

A non-zero fi s implies a non-zero pair breaking stress on ud and ds pairs. Once CFL 
pairing is disrupted, an additional stress arises due to the different electric charges of the 
u and d quark. There are two scales that govern the response to fi s , the mass tjik of the 
lightest strange Goldstone boson, the kaon, and the gap A for fermionic excitations. When 
H s is equal to m K the CFL phase undergoes a transition to a phase in which the CFL order 
parameter rotates in flavor space and a kaon condensate appears If fi s > A gapless 
ermion modes appear in the spectrum and the CFL pairing pattern starts to get disrupted 
6, 7]. We will refer to the phase with kaon condensation as the CFLK phase, and the gapless 
phases as the gCFL and gCFLK phase. 

Gapless fermion modes in weakly coupled pair condensates cause instabilities in current- 
current correlation functions j^l, 9] and these instabilities have been found in the gapless 
CFL phase llO , llll, . We have suggested that the instability is resolved by the formation 

nnn 

of a non-zero Goldstone boson current [13|, 114], |15j. In our previous work we analyzed the 
Goldstone boson current in the kaon condensed CFL phase, but we restricted the pairing 



ansatz to a flavor rotation of the CFL state and did not properly implement charge neu 



trality in the Goldstone boson current state Q, [jjjl. In 15|] we allowed the CFL pairing 
pattern to be disrupted and constructed a charge neutral solution, but we did not allow 
kaon condensation. 

In the present work we study Goldstone boson currents in a kaon condensed CFL phase. 
We use a very general ansatz for the current, allow the CFL pairing pattern to be disrupted, 
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and implement electric and color charge neutrality. We also compute the magnetic screening 
masses. Our work remains restricted to states in which the Goldstone boson current is small 
compared to the gap. This is expected to be the long as fj, s is close to the onset 

for gapless fermion modes in the spectrum. The opposite limit of a current which is large 
compared to the gap in analyzed in studies of the Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) 
state in three flavor QCD Q, 17, 18, Q, 



II. EFFECTIVE LAGRANGIAN AND DISPERSION LAWS 

Our starting; point is an effective lagrangian for gapped fermions interacting with gauge 
fields 



21 



22] 



C = Tr (xl(iv -d-fi L - A e Q) XL ) + Tr (x R {™ ■ d - fi R - A e Q) XR ) 
- iTi (xIxlXv ■ (d - iA T )X^) - *Tr ( Xr XrYv ■ (d - iA T )Y^) 

~ o A k (xLXLtkabtkij - XRXRtkabekij + fc.C.) . (1) 



2 a, b ,ij,k 



The notation agrees with our earlier work 15|] . x°lr are left/right handed fermions with 



color index a and flavor index i, A^ are SU(3)c color gauge fields, and fx L = MM' /(2/j,), 
fi R = M^Mj (2/i) are effective chemical potentials induced by the quark mass matrix M. 
The trace part of fi L ' R corresponds to a shift in the baryon chemical potential and will be 
neglected. The matrix Q = diag(|, — |, — |) is the quark charge matrix and A e is an electro- 
static potential. The fields X, Y determine the flavor orientation of the left and right handed 
gap terms and transform as X -> LXC T , Y -> RYC T under (L, R) e SU(3) L x 5^(3)^ 
and C G SU(3)c, and (k = 1,2,3) are the CFL gap parameters. From the lagrangian 
given in equ. ([1]) we can read off the left handed Nambu-Gor'kov propagator, 

where p = v ■ p — /x, with the Fermi velocity v. The right handed propagator is obtained 
by replacing (L) with (R). The components of the propagator are matrices in color-flavor 
space. We use a basis spanned by the Gell-Mann matrices X A (A = 1, . . . , 8) and A = y§ 1. 
In this basis 

A AB = -±A ab e ija e rsb \f r \f s , A ab = diag(A 1; A 2 , A 3 ) ab . (3) 



(2) 



The right and left handed gauge field vertices are (for Y = X^) 



X(l),v,ab = |Tr [X A {fi L + A e Q)X B + i\ A X B X v (d - iA T )X ] 
X(r),v,ab = |Tr \\ A (ji L + A e Q)X B + iX A X B X^ v {d - iA T )X 



(4) 
(5) 



Instead of using A3 and As it will turn out to be convenient to use the following linear 
combinations, 

A/ = A 3 H — t=A 8 , X11 = A 3 — V^Ag. 



V3 



We remark that [A/, A7] = 0. We assume a maximal kaon condensate 

(Z7T \ 



(6) 



(7) 



and make the following ansatze for A 0T and A T |23|, |24( , 

1 r 



A 



0T 



X\(i L + A e Q)X + X{fi L + A e Q)X^ + AjX! + A 7 X 7 



(8) 
(9) 



This particular flavor structure leads to diagonal matrices for XAj^X^ and X^A^X. For 
A 7 = the ansatz ([9]) for A T is equivalent to setting 



with 



U = exp 



-]■ x c 7 X H 



X = VU£ K oU^ 



V = exp 



-J- x c/Aj 



Co 
>/6 



A 



(10) 



(11) 



and taking the solution for A T given by A T = ^(X^dX + XdX^). For A 7 ^ the com- 
putation of dispersion relations using the ansatz (TTUT) is not straightforward because the 
x-dependence does not drop out. Therefore we will use the ansatz given by equs. (171) and 
(Q as our starting point. In the end we will see that the condition A 7 = is indeed fulfilled, 
which means that the ansatz for the gauge field (Q is equiva 
stone boson current fflQj) . We remark that the ansatz of Ref . 
c 7 7^ 0, while the ansatz of Ref. [3] corresponds to c , c 7 ^ 0, cj = 0. 

We find that the tadpole diagrams with an external ^4i, 2,4,5,6,77 vanish identically for any 
values of the parameters of our ansatz. Therefore the corresponding neutrality conditions 
are automatically satisfied for the system under consideration. 



ent to the ansatz for the Gold- 
corresponds to c = Cj = 0, 
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We find the following dispersion laws in the left handed sector, 

1*8 



H(L)2 



H(L)4 



"(i)5 



£(L)6 



9 -i 7 + (2 C/ + c )^ 



\ 



-^ + i 7 +(2c / + co)^ + 



C(L)s = ^ + i/-y-A e - (2(c 7 - c ) + 3c 7 ) i , ; 



2 

^P + T - T - T + (2c,-c 7 )^j +A 2 , 



3/x 



r -^/ + Y + ^-(2( C/ -c ) + 3c 7 )^ 



+ 



P + ^-^-^-(2c/-c 7 



A 7 A 



f • J 



f + i. + ^-A e -(2( C/ -co)-3c 7 )^ 



+ 



P- T - T + T + (2c, + c 7 )^ +A|, 



■^-^-^+^-(2(^-^-30,)^ 



+ 



\ 



' /x s A 7 A 7 u- A 

P- T - T + T -(2C + C 7 )^ + A 3 , 



e (i)7,8,9 = complicated. 



(12) 



The final three ("complicated") dispersion laws cannot be expressed in simple form. They 
are electrically neutral and do not have much impact on the dynamics, since they never 
become gapless for the values of fi s which we will consider. The dispersion laws for the right 
handed sector can be obtained from the left handed ones by substituting A 7 — > — A 7 and 
c 7 — > — c 7 . We note that for cj = c and c 7 = the current does not shift the energy of the 
electrically charged modes, which is similar to the case without kaon condensate. 



III. THERMODYNAMIC POTENTIAL 



In order to evaluate the thermodynamic potential and study the phase structure we will 
make a number of approximations. In QCD gluons are dynamical fields, and gluon exchanges 
determine the value of the gap parameter. In the following we will treat the gauge fields and 
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FIG. 1: Fermion gaps in the kaon condensed CFL phase as a function of the effective chemical 
potential fj, s . The figure shows A2/A(0) (upper curve) and Ai/A(0) (lower curve) as functions of 
/i s /A(0) for j = 0. Within our approximations A3 = A2. 



the gap as classical mean fields. The mean field partition function reproduces low energy 
QCD fermion propagators as well as n-point functions of the currents and gauge fields to 
leading order 0(a®) in the strong coupling constant and to 0(p%) in the effective chemical 
potential. 

In the mean field approximation the thermodynamic potential is given by 



(8c? + cl + 6c?)a*V f 



247T 2 



2tt 2 



1=1 



where the first term is the mean field potential for the gap field, the second is a tadpole 
term that arises from integrating out anti-fermions 251 ]. and the third term is the one- loop 
effective potential of fermions coupled to background gauge and gap fields. We will keep 
only the leading order terms in A//x in the thermodynamic potential. In particular, we have 
written the integral over the fermion dispersion relations in terms of the density of states, an 
integral dp over the momentum transverse to the Fermi surface, and an angular integral dv. 
We have included the A\ term which is formally of higher order, but necessary to determine 
a unique solution to the electrical neutrality condition. 

The integral over p is ultraviolet divergent and requires a cutoff A. This divergence can 
be absorbed into the coupling constant (7(A). In the following G refers to G(A = p). The 
constant G fixes the magnitude of the gap in the chiral limit, which we denote by A(0). 
We shall study the thermodynamic potential as a function of the dimensionless variables 
/i s /A(0) and j/A(0). 
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The integrals over ei_6 are easily evaluated. The result is 



ft = —(A? + A2 + A3) + ^ 1 - 7;p J + ft 



12tp 



+ 



+ 



24vr 2 

/ 2i/ n 8 ~ jc 7 j\ 

{ Al >^r>Y- A7 >^ {2ci+Co) 6) 

t? ( a ^ . Ai . A7 3/i s ~ A 7 j 
(2( C/ -c ) + 3c 7 )^ 

+ F ( a 3 , ^ - ^ + il, ^ + I, + - A e , (2c 7 + c 7 ) 3 -. 
V 4 3 2 4 2 v v -r 7 ; 4 , 



(2( C/ -c )-3c 7 )^J+(A 7 



where ft is the integral of the complicated modes, 

,2 , 9 



-A 7 ,C 7 



-or) 



ft 



2tt 2 



The function F(A, a\, 02^1, J2) is defined by 



F(A,ai,a 2 , ji,j 2 ) 



127T 2 



2/i 



3A 2 (l + 21og(^))+6a 2 + 2j 2 



|e(a 2 - j 2 - A)- [A(j 2 ) (A(j 2 ) 2 + 3A 2 ) 



24vr 2 



+3A 2 (a 2 -j 2 )log 



32 

'0-2-32- KJ2) 
a 2 ~ 32 + K32) 



+ (J2 



-32, 



(14) 



(15) 



(16) 



with A(j 2 ) = y (a 2 — j 2 ) 2 — A 2 . Here, a x and ji refer to the shift in the minimum of the 
dispersion relation due to the gauge potentials and the current, and a 2 and j 2 denote the 
shift in the energy. The equations for color and electrical neutrality and the equation 
dn/dA 2 = dtt/dA 3 = are solved by 



Aj = A 7 = 0, 



A 3 = A, 



^3 — ^2, 

is 



A, = max ( 0, -A 2 + + (2|c 7 - c | + 3|c 7 |) ^ ) . 



(17) 



which greatly simplifies the expression for the thermodynamic potential. The remaining 
equations for Ai, A 2 and j are solved by numerical minimization of the thermodynamic 
potential. 

First let us consider the case j = 0. Fig. [1] shows the solution of the gap equations. 
We notice that even for small fi s there is a small splitting between Ai and A 2 3 . The 
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FIG. 2: Thermodynamic potential of the kaon condensed CFL phase fi/(^ 2 A(0) 2 ) relative to the 
thermodynamic potential of neutral non-interacting quark matter as a function of /i s /A(0) for 
j = 0. For comparison we also show the thermodynamic potential of the gCFL phase using the 
same approximations (dotted line). 

thermodynamic potential is shown in Fig. El where we also compare to the result for the 
CFL-gCFL phase without a kaon condensate. We find that €4 becomes quadratic in p for 
H s > = 1.35A(0), and that the mode e 2 becomes gapless for fi s > /4 2 ) = 1.607A(0). In 
contrast to the CFL-gCFL phase transition, the phase transition at fj,y) is first order. At 
H s = there is a substantial drop in A x and a small increase in A 2 3 . The thermodynamic 
potential is smooth at fffl but has a kink at fi^\ 

We observe that near fi^ the thermodynamic potential of the gCFLK phase is very 
close to the thermodynamic potential of the gCFL phase (fl g cFLK < QgCFL by a very small 
amount). It is not entirely clear why this is the case. Near jj,^ the gauge potentials become 
quite large and our approximations are breaking down. Also, since some of the gaps become 
quite small it is not clear whether the correct CFLK state is a simple flavor rotation of the 
CFL state. We note that near jj,^ the thermodynamic potential in the gCFL phase plotted 
in Fig. [1] is lower than the result of Alford et al. [7j . This difference is also due to higher 
order terms. In the gCFL phase the thermodynamic potential in equ. (fT3|) differs from the 
functional used by Alford et al. by terms of 0(/4). 

Next we study the effect of a Goldstone boson current. We first consider the region fi s ~ 
fi^p, which was also investigated in [k], [yj]. If we take into account the condition of electrical 
neutrality we find that the expectation value of the current becomes extremely small. The 
competition between the A\ and fi 2 j 2 terms leads to a minimum of Q at j ~ (|/i s — A 2 ) 3 / /i 2 , 
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. 0008 
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FIG. 3: Thermodynamic potential as a function of the Goldstone boson current near /4 • The 
figure shows = 0)]/(/i 2 A(0) 2 ) at y. s = 1.600A(0) (dotted), /i s = 1.605A(0) (continuous) 

and Li s = 1.610A(0) (dashed). 

corresponding to a contribution to Cl of the order (|// s — A 2 ) 6 //z 2 . A detailed calculation of 
the current is presented in the appendix. Since the current is so small, we shall neglect its 
effect in the remainder of this section. 

Now consider the region /x s ~ fi™. Here the optimal value for c-j turns out to be zero, 
because c 7 does not contribute to the shift of the energy of e 2 . If we set cj = c the current 
shifts only the energy of the electrically neutral mode. This should be the preferred current 
if electrical neutrality is enforced. We set cj = 1.2 and Cq = 1, which leads to an even 
stronger minimum of the thermodynamic potential at finite j. 

Fig. [3] shows the thermodynamic potential as a function of j for different values of fi s in 
the vicinity of fi s = 1.6A(0). The kinks come from a discontinuity in A 12 (j), similar to the 
discontinuity in Ai i2 (/x s ) shown in Fig. [TJ We observe that the thermodynamic potential 
develops a nontrivial minimum if fj, s > 1.605A(0), which is slightly below the value of /j, s 
where the mode e 2 becomes gapless at zero current. We remark that the window where the 
thermodynamic potential has a nontrivial minimum is rather small. For fi s slightly above 
1.615A(0) we find that the global minimum is again at j = 0. We note, however, that the 
some of the gaps are very small, the current at the non-trivial minimum is rather large, 
and our approximations are not reliable. Since the currents are large the correct ground 
state may well be of the LOFF type and involve multiple currents and complicated crystal 
structures as in 26j. In addition to that, inhomogeneities of the amplitude of the order 



parameter may play a role 27, 
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m Af/( m MU.=o) 



1 




FIG. 4: Magnetic gluon screening (Meissner) masses in the kaon condensed CFL phase as a function 
of the effective chemical potential fi s . The figure shows squared Meissner masses rn 2 MU = m 2 M22 = 
m M44 = m M55 (continuous), m 2 Mm = m 2 M ^ (dashed), m 2 M77 (dotted), rn 2 M ( 2 ) (dash-dotted) in 
units of m^|^ s=0 = (21 - 8 log 2)g 2 /j 2 /(W8it 2 ) as functions of // S /A(0). 

IV. MEISSNER MASSES 



In this section we shall study the stability of the CFLK and gCFLK phase with and 
without a Goldstone boson current with regard to chromomagnetic instabilities. At leading 
order in a s the screening masses are given by a one-loop diagram involving Nambu-Gor'kov 
propagators. We have 

(™mU - 27f2 <W + lg7r2 1™ hm y d P y ^ y ^ 

x Tr [G+{p)Vt L)a G+{p + fc)Vj )6 + <T (p)Vfo a G-(p + k)V> Lyb 

+ S + (p)\/ ( l )a S-(p + k)V( L)b + S-(p)V7 L)a S + (p + fc)Vj )6 + (i2 <-> L)} ,(18) 

where we have denned the vertices 

(Vj )o )^ = ^r[\ A \ B X\ T a X^]v\ (VfojAB = -lTr[\ B \ A X\ T a X^\ 

{Vl R)a ) AB = \Tv[\ A \ B X^\lX}v\ {V( R)a )AB = -\Tv[\ B \ A X^\ T a X}v\ 

with v l = —v l . We remark that right handed and left handed quark propagators are identical 
because of A 7 = c 7 = 0. 

We find that the mass matrix in the 3-8-sector is not diagonal, the mixing angle being 
equal to ir/6. We denote the corresponding eigenvalues by m 2 M ^ and m 2 M ^ 2 y We find that 
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2 

m Mil 


2 

m M66 


2 

m M77 


2 

m M{2) 


transv. 


0.16 


0.59 


1.09 


0.21 


long. 


0.31 


0.87 


1.05 


0.46 



TABLE I: Meissner masses at the minimum of f£(j) for /i s = 1.605A(0), in units of m| / | /Js= o. 
some of the Meissner masses are equal, 

2_2_2_2 2_2 (1Q\ 

m Mll — m M22 — m M44 _ m M55> m M66 — m M(l)- \ 1J ) 

Fig. H] shows the Meissner masses as functions of /i s for /i s < //( 2 ) = 1.607A(0). The kink in 
mMii at [x s = fj/f' results from the kink in A e . We observe that the Meissner masses in the 
gapless CFLK phase are real, even in the absence of a current. 

For /i s > /4 2 ) and zero current the Meissner masses m Mll = m\ l22 = m 2 MM = m M55 as 
well as m 2 M ^ receive a negative contributions proportional to — — 4Af. This is the 

usual chromomagnetic instability for zero current. In the presence of a finite current J we 
may decompose the Meissner masses into a longitudinal and a transverse component, 

{m 2 M f = m M± (S ij ~ ff) + m 2 M] 0. (20) 

The values of the Meissner masses at finite current just above the phase transition at fi s = 
1.605A(0) are shown in Tabled! We note that all these values are positive. We note that for 
fi s > 1.615A(0) there is no Goldstone boson current and the magnetic screening masses in 
the CFLK phase are tachyonic. This is an unusual situation because stability with respect 
to a Goldstone boson current implies that the phase is stable with respect to electrically 
neutral fluctuations of the gauge field. We will comment on this regime in the conclusions. 

V. CONCLUSIONS 

We have studied the stability of a kaon condensed CFL phase with respect to the for- 
mation of a Goldstone boson current. We have computed the thermodynamic potential and 
the gluon screening masses as a function of the effective chemical potential /i s = m 2 s /(2pF)- 
Our starting point is an effective lagrangian of fermions coupled to a mean-field gap term 
and background gauge potentials. Our methods are reliable provided the gauge potentials 
and currents are smaller than the gaps. There are three important scales that appear in 
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the problem, /x s = m K is the onset of kaon condensation, /i s = /iW = 1.35A(0) is the point 
where a charged fermion mode becomes (almost) gapless, and fi s = ^ = 1.61A(0) is the 
point where a neutral fermion mode becomes gapless. The CFLK-gCFLK transition at fi^' 
is smooth while the transition at /j,^ is first order. 

We find that the kaon condensed CFL pha se is more stable with respect to chromo- 
magnetic instabilities than the CFL phase ^{J . 0|. The gluon screening masses are real 
in the regime fi s < /4 2 \ despite the presence of an almost gapless mode. There is a weak 
Goldstone boson current instability for fx s > iJp . The magnitude of the current is suppressed 
by the constraint of electric charge neutrality. 

We find that the gluon screening masses become imaginary in the CFLK phase for fi s > 
fj,@\ There is a Goldstone boson current instability that develops at this point, and the 
screening masses in the presence of a Goldstone boson current are positive. However, the 
Goldstone boson current phase only exists in a very small window above /4 ■ Beyond 
that point the CFLK phase is stable with respect to the formation of Goldstone boson 
currents, but the gluon screening masses are imaginary. This is somewhat puzzling since 
the Goldstone boson currents are equivalent to non-zero gauge fields, see equ. @. The 
difference between the second derivative of the thermodynamic potential with respect to 
the Goldstone boson currents and the screening masses is that the former are computed 
at constant (zero) electric charge, while the latter are computed at constant electro-static 
potential. In this sense stability with respect to Goldstone boson currents is the physically 
relevant criterion. The fact that the Meissner masses are imaginary may nevertheless imply 
that other instabilities are present. In this regime some of the gaps are small compared the 
background fields and it may be necessary to consider states with multiple currents, such as 



the crystalline LOFF state [18|, |2Q, [26| . 
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APPENDIX A: CHARGE NEUTRALITY AND GOLDSTONE BOSON CUR- 
RENT FOR n s > 

In this appendix we study in more detail the current instability in the regime fx s > jj^ ~ 
4A/3. For simplicity we only consider a pure hypercharge current with V = 1. Including 
the baryon current will only change some numerical coefficients. Our analysis extends the 
analysis of Q] to properly implement charge neutrality in the Goldstone boson current 
state. 

The instability is driven by the lowest fermion mode in the spectrum. The dispersion 
relation is given by 

A (p-lo) 2 3 U - , A ,x 

UJ p = Jt ~ 2A ~ 4^ + 6 ~ l V ' 3Kl ( ^ 

where A e is the electron chemical potential, fi s = m 2 s /{2p F ) and l = (/^s — v • ]k)/4. The 

contribution of gapless modes to the thermodynamic potential is 



^ = j dp ltr (A2) 



where the factor 2 is a degeneracy factor and dv is an integral over the Fermi surface. The 
charge density of gapless modes is 

The electron contribution to the thermodynamic potential and the charge density is 



A 4 A 3 

IZ7T 2 37r / 



We introduce dimensionless variables 

X "aA' aA ' K " aA ' (A5) 

where a is a numerical coefficient that depends on the parameters of the effective theory and 
the number of currents that are turned on. For a pure kaon current 

' (A6) 



15 2 4<' 



where c n = (21 — 81og(2))/36 is the numerical constant that appears in the weak coupling 
result for the pion decay constant f n and v\ = 1/3 is the square of the Goldstone boson 
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velocity. The thermodynamic potential and charge density can be written as 

/i 2 A 2 



Q(h, h e , x) 



7l z 



(cf h _ he (x)-cx), 



/i 2 A 

q(h,h e ,x) = — ^- (Kg h ^ he (x) - K e hfj 
f h {x) = x 2 - - \(h + x) 5/2 Q{h + x )-(h- xf /2 Q{h - x) 

T. L 



with 



9h ( x ) = -\(h + x) 3/2 Q(h + x )-(h- xf /2 Q{h - x)] . 

X 1 J 



The numerical coefficients are given by 

2 fi 2 A 2 



and 



C 



a 



— (-) 

3 • 4 5 I ^ ) ' 



K 



K P 



IOC 



i6a 



(A7) 
(A8) 

(A9) 
(A10) 

(All) 
(A12) 



We first solve the neutrality condition in order to find h e = h e (h,x). We have h e — if 
h + x < 0. For h + x > we can use K e <C 1 in order to write 



with 5/i e <^ /i e . We get 



h e = h + x + 5h e 



^) x 2 '\h + x) 2 . 



(A13) 



(A14) 



Inserting this result into the energy density gives 



/x 2 A 2 



n 



c 



x 



-Sh 5 J 2 e(Sh e )] - CJh + x) 4 + ...}. (A15) 

X J J 

The contribution from gapless fermions is proportional to x 2 ^ 3 for small x so the instability 
is still present. If K e , C e <C 1 the minimum of the thermodynamic potential is determined 



by the balance between the electron term and the current contribution. We get 



2C e;3 /A 2 \ 3 



C 



(A16) 



This implies that a Goldstone current is formed, but the magnitude of the current is sup- 
pressed as compared to the result given in [l3|, [14| . 
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